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The non-linear current-voltage characteristic of a tunnel junction between two Luttinger systems 
is calculated for an interaction with finite range. Coulomb blockade features are found. The dissipa- 
tive resistance, the capacitance and the external impedance, which were introduced ad hoc in earlier 
theories, are obtained in terms of the electron-electron interaction. The frequency dependence of 
the impedance is given by the excitation spectrum of the electrons. 

PACS numbers: 73.40.Gk, 72.10.Bg, 72.10.-d 
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The Coulomb blockade effect in the non-linear current 
voltage (I-U) characteristics of mesoscopic tunnel junc- 
tions has been the subject of many theoretical and 
experimental investigations ||Q during the past decade. 
Basically, due to the repulsion between the electrons, 
tunneling is suppressed for voltages below Uc = e/2C, 
and temperatures smaller than Tc = Ec/ks {ks Boltz- 
mann constant, e elementary charge, C capacitance). 
The quantity Ec = ellc is the charging energy. 

In the semi-phenomenological theory of the phenome- 
non 1^ the tunnel junction is modelled by a capacitance 
and a tunnel resistance Rt- An impedance Z(lu) is in- 
cluded into the circuit. It represents the coupling of the 
tunneling particles to a reservoir of Bosonic degrees of 
freedom. They guarantee incoherence between different 
tunneling processes. When Z{0) = R — 0, the current 
voltage characteristic is linear, I{U) = U /Rt- For R ^ 0, 
I(U) cx [/2«/«K+i .^^hen U <s:Uc {Rk von Klitzing con- 
stant), and I{U) w {U~Uc)/Rt when Uc <.U oo. It 
must be emphasized that the shift of the linear behavior 
of I{U) by Uc is the important characteristic feature of 
the Coulomb blockade for i? ^ oo. 

In this paper, we present a microscopic theory of the 
effect for a one-dimensional (ID) tunnel junction. The 
parameters introduced in the above mentioned theory by 
ad hoc assumptions are deduced consistently, and in a 
natural way, from the interaction between the electrons. 

Two semi-infinite (ID) systems of interacting electrons 
described within the Luttinger approximation Q are cou- 
pled by a tunnel junction. The interaction potential be- 
tween the electrons is assumed to have a finite, non-zero 
range. The tunneling current as a function of a voltage 
applied across the junction is obtained. 

The charging energy is found to be the interaction po- 



tential at zero distance, and the dissipative resistance is 
given by the spatial average of the interaction potential. 
The spectrum of the elementary excitations determines 
the impedance of the circuit. In order to explain the lat- 
ter no additional 'environmental modes' are needed. The 
crucial point is that the above mentioned asymptotic be- 
havior of /([/) for large U appears to be directly related 
to the finite, non-zero range of the interaction. The latter 
implies that the dispersion relation of the elementary ex- 
citations of the electron system becomes u:{k) « VF\k\ in 
the short wavelength limit [Q. Our results show that in 
ID a charging energy, and, in turn, a capacitance can be 
defined in the Luttinger model provided the interaction 
has non-zero range. 

We consider the Hamiltonian H = Ho+Ht+Hjj. Here, 
Hq = nj^p consists of the Hamiltonians of the two 

disconnected electron systems, which extend from —L to 
and from to L (L ^ oo), respectively. The tunnel 
junction (at a; = 0) is described by Ht, and Hu is the 
energy contributed by the external voltage. 

The electrons are described by the Bosonic Luttinger 
Hamiltonian S, (j = 1,2) 
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Boundary conditions are assumed such that the origi- 
nal Fermion fields vanish at x = 0, ±i. This introduces 
additional (quadratic) off-diagonal terms in the Hamilto- 
nian . We neglect them here for the sake of simplicity, 
since they do not affect the final results qualitatively. 
The dispersion law uj{q) = qv{q) with the wave num- 
ber dependent velocity v{q) = vp[l + V{q)/TTVF]^^^ , re- 
fiects the Fourier transformed of the interaction potential 
V{q). For the latter, we assume a 3D screened Coulomb 
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potential with the range projected onto a quantum 
wire of diameter d. Depending on whether d 
or a^^ ^ d the interaction is exponentially (Luttinger 
limit) or algebraically (cx x~^, Coulomb limit) decaying, 
respectively. For small q (= mr/L, n integer) we have the 
charge-sound excitations with the renormalized velocity 
v(0) = vp/g characteristic of the Luttinger system, with 
= [1 + V{0)/ttvf]^^^- For large q, we find the exci- 
tation spectrum of the non-interacting electrons, due to 
the finite range of the potential. The velocity associated 
to AA^ (see below) is vn = v{0)/g. 

The operators j'^^^^ ^ 7^^' are related via a Bogolubov 
transformation to the Fourier components of the 
phase fields, bq^^ = coah{ipq)jq^'' — sh{ipq)jq^^\ 



They define (right moving) Fcrmion fields 



1 g-ii3<,^''gi7ra:AAf<^V-f-g«*'^'(2:) 
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with the density operators p-^ {x) = ^'^-'^(x)*^^(a;). 
The functions Lpq contain the above dispersion relation, 

exp {-2ipq) = uj{q)/vFq- 

(i) 

The variables i?q , defined modulo 27r, are conjugate 
to the number operators, ■dQ\ AN^^^ = i, 



(4) 



The latter represent extra electrons in the systems on the 
left and the right hand sides of the junction. 

The above boundary conditions imply that the corre- 
sponding left- and right moving parts are not indepen- 



dent but ^-^^(0;) = ~^l>{~x), ^^^>{x + 2L) = ^^^\x), 
and either one of the two alone suffices to describe the 
system. The tunnel Hamiltonian in terms of the latter is 

00 
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(0)*g'(0) 
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By inserting the above Bosonized form (|^) one obtains 



Ht = H+ + H^r, with 
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iJf ^ = exp <( ±* ^ (-f )■'" - ci>(^) (0)] \ . (6) 

The electrostatic energy of the external voltage that 
is assumed to drop only at the tunnel junction, U{x) = 
U [e(x) - e(-a;)] /2 (e(x) Heavyside function), is 



Hu^-e I ^ U{x)p{x), 



(7) 



with p{x) = p(i)(x)e(-a;) -|- p'-'^'^ {x)Q{x). In the dc 
limit, it has been shown that the current- voltage rela- 
tion is independent of how the voltage drops | p^ , p^ . 
Only the voltage drop between x —00 and x —> 00, 
which is assumed to be fixed by an external "battery", 
is important. Therefore, interaction induced rearrange- 
ment of charges in the presence of the impurity is unim- 
portant for the present calculation. By inserting the 
above relation (^ between p^"* and AiV^^', noting that 
p(^")(x) = p^^\x)+p^i\x) = p'^^\x)+p^^\-x), and using 
{L) = (-i) one obtains 
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The current operator is defined as in previous works 
/ = ie[H^ ~ H^] jl4j. Fermi's golden rule with Ht as a 
perturbation fin] yields the average current 
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-^^^] / dte'^^^^e-'^^^*) (9) 



{f3 inverse temperature). This equation is independent 
of the boundary conditions applied. The thermal equi- 
librium correlation function 

(t)^ J2 { [*^'^ " (0)) (10) 



Wg(t) = 



is evaluated with respect to Hq- The ^^^'^ {t) evolve in the 
interaction picture with respect to Hq + Hjj. This gives 



W,it) = I du^^ X 



(I — cos ([x)t))coth 



I3lo 



- i sin (cut) 



(11) 



The key quantity is the spectral function J{lo) which is 
directly given in terms of the dispersion law of the charge 
excitations of the system 



J{u;) 
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where q{uj) is the inverse of the dispersion function Lo{q). 
Equation (^2|) establishes the main result of this work in 
the sense that the quantity on the left hand side, which 
was already present in the former theory is now di- 
rectly related to spectrum of the elementary excitations 
of the electron system Q . 

The impedance function Z{lu) = J(ui)/uj — 2 for the 
interaction potential of the Luttinger limit is shown in 
Figure . Qualitatively similar curves are obtained for the 
Coulomb limit. For w ^ 0, Z{0) = 2(5-1 - 1) = R/Rk- 
For large lo, Z{uj) tends to zero, i. e. J{(^) ^ 2w, the 
limit of non-interacting electrons. The peak in Z(uj), 
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which appears for smaller values of the interaction pa- 
rameter g, is a consequence of the inflection point in the 
spectrum of the elementary excitations at LUp (the fre- 
quency of the maximum in q'{uj) 0). When the screen- 
ing length decreases, LUp, and consequently the frequency 
beyond which Z(u)) vanishes, increase proportional to a 
(Fig. ), for a ^ oo, Z{u;) = 2{g~^ - 1). 

In order to calculate the current voltage characteristic 
we have to Fourier transform exp [— Wg(t)]. This is done 
by using the relation between Wi{t), and the Fermi dis- 
tribution f{E) 
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The result is 
/([/) 



eA^ 
27r 
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/OO 
dEf{E)[l- f{E + eU)] 
-OO 



^D{E)D{E + eU) 
with the tunneling density of states 

/OO 
dt COS (M)e-^«(*)/^ 
-OO 



(14) 



(15) 



It is nothing but the local density of states at the position 
of the tunnel junction |M,n]],nffl. 
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FIG. 1. The impedance function J{ui)/uj — 2 of the Lut- 
tinger limit for different interaction parameters g. 

The connection with the earlier results j|] is established 
by observing that the I{U) can also be written as 

-l3eU 
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dE 
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dE'f{E)[l-f{E')] 
xP{E + eU - E'), (16) 



with the probability density for a bulk excitation of en- 
ergy E, 
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^_^^,Et^-[W,{t)~W,(t)]_ 



(17) 



The tunnel resistance is Rt = 2a;^/e^A^7r [T^ . 

In P{E) the role of the 'electromagnetic environment' 
is now played by the excitations of the interacting elec- 
trons. The result for zero temperature. 



^(^) = ^ / dE{eU - E)P{E), 
eRt J a 



(18) 



is shown in Fig. for different interaction strengths. 

For smah voltages I{U) « {1/ Rt)U'^/3-^ . By com- 
parison with the corresponding limit of reference |^ , we 
recover the dissipative resistance in terms of the interac- 
tion 1I0|,|T1, 



R 
2R^ 



1 

- - 1. 
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(19) 



For U much larger than the range of P{E) eq. ( p^ ) 
becomes 



I{U) = 



u- 



Ec 



(20) 



from which the charging energy is found by using our 
microscopic expression ( |l2|) for J{uj) 



En - 



dujZ(uj) = 2V{x = 0). 



(21) 



Taking into account the above mentioned off-diagonal 
terms induced by the boundary conditions docs not chan- 
ge the dependence of Ec on the interaction potential. 
Only the prefactor is altered 0|. 
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FIG. 2. The current voltage characteristic of a tunnel junc- 
tion between ID interacting electron systems for different 
strengths of the interaction, g. 

The capacitance found from Ec, C — 6^/41^(0), scales 
in the same way as the one, Cq, found previously by 
considering the frequency dependent conductance of a 
pure quantum wire of interacting electrons [Q. This 
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was obtained by a completely different philosophy start- 
ing from the classical theory of antennas. In fact, in 
the Coulomb limit, C/Cq = 0(1). Following this ap- 
proach, one could consider the present problem as being 
the quantum equivalent of two only very weakly coupled 
wires. From the comparison of their AC properties with 
those of the quantum system considered here, one would 
get the same scaling law for the capacitance. 

The capacitance responsible for the Coulomb blockade 
is naturally given here by the total capacitance of the 
circuit, i. e. the two wires and the junction. Apparently, 
in our ideal ID model, the junction as such does not con- 
tribute significantly. It is only needed for detecting the 
Coulomb blockade induced by the latter. In the Coulomb 
limit, C = neeod w O.OSfF, for a quantum wire fabri- 
cated in a semiconductor heterostructure, AlGaAs/GaAs 
(e « 10, d ~ lOOnm), for instance, with only one subband 
occupied. The charging energy is Ec ~ 2.3meV, which 
corresponds to a critical temperature of Tc = 30K. By 
decreasing the width of the wire, which seems achievable 
with present days' technology, one should be able to in- 
crease Tc close to room temperature. Thus, a one-mode 
narrow quantum channel with a weak link appears to 
be the 'ultimate device' for observing Coulomb blockade 
effects 0. 

From the AC properties |^ , we find also an inductance 
L of the electron system, due to the presence of the reso- 
nance in the AC conductance at LOp for strong interaction. 
It is determined here by the ratio Ec /cOp = y^L/C, and 
is reflected by the resonant behavior of the impedance 
(Fig. ). This is, however, of negligible influence on 
the behavior of I{U). Since the integral weight of the res- 
onance is small compared with the charging energy, the 
total integral over Z{uj), and also uip <C Ec {Ec ~ bOOup 
for the above single mode wire), the steps in the deriva- 
tive of the DC current voltage characteristic predicted 
earlier as a signature of a (5-function like resonance in 
Z(uj) cannot be obtained by using the present model. 
This prevents the inductance to be detected directly here. 

Our above results remain qualitatively true for a non- 
linear dispersion relation, beyond the Luttinger model, 
since the influence of the interaction is negligibly small 
for wave numbers much larger than the inverse of the 
range of the interaction. This is indicated by the large q 
behavior of the dispersion with interaction in the random 
phase approximation |p^ . 

For many channels occupied some of the above conclu- 
sions are modifled. Matveev and Glazman ||l^ treated 
the tunneling for a quasi-lD quantum wire with many 
(TV) channels, but for zero range interaction. They flnd a 
crossover in the zero frequency impedance from eq. ( |l9| ) 
to the non- interacting limit (i? — s- 0) when N oo. 
In that model, one cannot obtain the asymptotic linear 
behavior of I{U), even for finite N. Therefore, we con- 
clude that our above result - that using the Bosoniza- 



tion method the Coulomb blockade is closely related to 
the finite range of the interaction - remains also valid for 
many channels, though the quantitative behavior of the 
charging energy, and thus the capacitance, will be prob- 
ably changed. We expect a contribution of the tunnel 
junction to the capacitance proportional to the channel 
number. This will be the subject of future studies 

In summary, we obtained the non-linear current volt- 
age characteristic for a model of two ID quantum wires 
of interacting electrons connected via a tunnel junction. 
The features of the Coulomb blockade phenomenon were 
found. The charging energy, the capacitance, and the 
inductance of the circuit were given in terms of the inter- 
action potential. In order to obtain the Coulomb blockade 
in the Luttinger liquid model, it is necessary to assume 
that the interaction is of finite, non- zero range. The role 
of the modes of the environment is played by the elemen- 
tary excitations of the interacting electrons. 
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